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The electronic structure of crystalline CdTe, CdO, Q-Te02, CdTeC>3 and CdaTeOs is studied by 
means of first principles calculations. The band structure, total and partial density of states, and 
charge densities are presented. For cv-Te02 and CdTeC>3, Density Functional Theory within the Lo- 
cal Density Approximation (LDA) correctly describes the insulating character of these compounds. 
In the first four compounds, LDA underestimates the optical bandgap by roughly 1 eV. Based on 
this trend, we predict an optical bandgap of 1.7 eV for CdaTeOe- This material shows an isolated 
conduction band with a low effective mass, thus explaining its semiconducting character observed 
recently. In all these oxides, the top valence bands are formed mainly from the O 2p electrons. On 
the other hand, the binding energy of the Cd 4d band, relative to the valence band maximum, in 
the ternary compounds is smaller than in CdTe and CdO. 
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INTRODUCTION 



Cadmium Telluride, CdTe, has been considered a pro- 
totype of II- VI semiconducting compounds for more than 
30 years. It has found important applications in 7-ray 
detectors, infrared windows, solar cells and other opto- 
electronic devicesi Having being the object of many of 
the early studies in semiconductor science, the difficul- 
ties found in obtaining high quality crystals hampered 
its research and applications in cutting edge electronic 
technology. Recently CdTe has received renewed interest 
due to the search for cheap technologies for mass pro- 
duction of solar cells that do not require high quality 
monocrystals. The fact that the fundamental absorption 
gap of CdTe lies in the region of maximum intensity of 
solar radiation makes it an important material for solar 
energy conversion on the earth. The solar cells based on 
heterojunctions of CdTe/CdS are a valuable option. 

The development of a complete technology requires 
complementary materials that could be used for sub- 
strates, coats, contacts, gates, etc. Among these ma- 
terials are the amorphous CdTe oxides obtained by r.f. 
sputtering, which have been extensively studied by opti- 
cal spectroscopy, X-ray photoemission spectroscopy and 
Auger electron spectroscopyi 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i 1:L i 12 i 13 These 
alloys are transparent insulators whose optical gap can 
be tuned according to the content of oxygen between 1.5 
and 4 eV, and can play a role in CdTe technology similar 
to that of Si02 in Si technology. Native oxides have also 
been identified on CdTe surfaces ^i^iLiLi^i However, 



given the difficulties produced by its amorphous struc- 
ture, the structural and electronic properties of these ma- 
terials arc mostly unknown. Even the nature of oxides in 
CdTe surfaces is still a subject of controversy^ 

The development of a thorough understanding of the 
electronic properties of the amorphous CdTe oxides re- 
quire a careful study of the crystalline compounds of Cd, 
Te, and O. These compounds provide information about 
short range order and local electronic properties that 
are determinant in the amorphous regime. The Cd-Te- 
O compounds can also provide valuable information to 
help understand the processes of oxygen interaction with 
crystalline CdTe, including adsorption and diffusion. 

Modern ab-initio simulation provides a powerful tool 
to study complex compounds and to predict struc- 
tural properties with reasonable confidence. Although 
there are many reports of ab-initio calculations of Cd 
and Te compounds, including c-CdTe, 20 i 21 i 22 i 23 l-CdTe^4 
c-CdQ 22 i 25 . 26 i 27 i 28 i 29 and TeO„ clustera 30 ! 31 . 32 , 33 , 34 :^, 
the crystalline Te02 and the ternary compounds remain 
largely unexplored. Only one ab-initio study of a Cd- 
Te-O system, which addresses the O adsorption on CdTe 
surfaces, is known to the authors* 1 - 9 - Despite being a prof- 
itable material for electromechanical and acoustooptical 
devices^ in the a phase (paratellurite) even crystalline 
Te02 has never been the object of a full electronic struc- 
ture calculation. 

The main difficulty for semi-empirical or ab initio simu- 
lation of Cd-Te-O compounds has been the large number 
of atoms per unit cell in the known phases j£L2iiSL2il (e.g. 
12 atoms in a-Te0 2 , 24 in /?-Te0 2 , 20 in Cd 3 Te0 6 , 40 



in CdTe0 3 , 224 in CdTe 2 5 ). Moreover Cd atoms have 
ten electrons in level 4d which lie in the range of the va- 
lence band of Cd compounds and increase significantly 
the computer resources necessary for the calculations. In 
certain cases, the effect of the Cd 4d electrons can be rep- 
resented with the help of pseudopotentials plus a partial 
core correction scheme; this approach has been successful 
for molecular dynamic simulations. 24 However, every use 
of this scheme needs to be tested, as the Cd 4d electrons 
are known to influence several important properties, e.g. 
the structure of the valence band, bandgaps, cohesive en- 
ergy, and lattice parameter iSl 

In this paper we report the first ever calculation of 
the electronic properties of crystalline a-Te02 (paratel- 
lurite), monoclinic CdTeC>3 (Cadmium oxotellurate(IV)), 
and CdaTeOe (Cadmium hexaoxotellurate). In order to 
analyze the nature of the chemical bond in these materi- 
als, and to evaluate the systematic errors inherent to the 
theory, we also analyze the electronic structure of the 
simpler semiconducting compounds CdTe and CdO. 

II. STRUCTURAL PROPERTIES 




FIG. 1: Structural unit of a-TeC>2. Generated with 
XCrysSDenii* 4 ^ 

Te02 is a mixed covalent-ionic insulator, the most 
abundant phase of which is a-Te02 or paratelluritc. 
Paratellurite structure has a tetragonal lattice with four 
formula units per unit cell. The space group has been 
identified as P4i2i2 and P43222. 4 4 In this work we have 
studied the electronic properties for the most recently 
reported^ At the short-range level, the structural units 
of paratellurite are Te04 trigonal bipyramids (tbp) with 
Te atoms at apex (see Fig. ^| . The corners of the tbps 
are occupied by two equatorial oxygen atoms (O eq ), with 
Te-O e<? bond length 1.88 A and angle O e «-Te-O e<? of 103°, 
and two axial oxygens (O ax ), with Te-O ax bond length 
2.12 A and angle O ax -Te-O ax of 168°. The tbps form 
a three dimensional network, interconnected through an 
equatorial and an axial oxygen. This short-range order 
can be understood as produced by Te sp 3 hybridization. 
This configuration consists of the sharing of four of six 
valence electrons, and the repulsion of four bonds by the 
lone pair of Te electrons^ Alternatively, the basic unit 
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FIG. 2: Atomic environments in CdTeOs. Generated with 
XCrysSDerM. 

can be considered as a TeC>6 distorted octahedron (see 
Fig. |SJ), which is obtained by adding to a tbp the next 
nearest pair of O at 2.86 A from the central Te atom. 
From this point of view the crystal structure of paratellu- 
rite is considered a deformation of the structure of rutile. 
This approximation has been used in previous bandstruc- 
ture calculations^ 5 *^ Ab initio calculations of TeC>4 and 
Te0 6 clustera 32 i 33 i 35 show that the Te-O ax bond is less 
covalent than Te-O e? , and the bond with the third pair 
of oxygen atoms is rather ionic. Recent studies^ii show 
that the dynamic properties of paratellurite are better de- 
scribed if the basic unit is considered as molecular TeO^ 9 , 
while Te-O ax bonds are understood to be intermolecular 
contacts. These bonds are important for the non linear 
dielectric properties 

Crystalline CdTeOa has been obtained in a cubic phase 
and several monoclinic phases i 3 ^^ The cubic phase has 
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FIG. 3: Atomic environments in CdaTeOg. Generated with 
XCrysSDenii*^. 



also been reported in oxidized CdTei^^i^^ Only one of 
the monoclinic structures has been totally determined^ 
This structure belongs to the P2i/c monoclinic crys- 
tallographic class, with the atoms occupying 10 non- 
equivalent positions. Figure |2] shows the first coordina- 
tion shell of all the non-equivalent atoms. The struc- 
tural units are distorted TeC>3 trigonal pyramids that are 
linked to each other with Cd atoms. This arrangement 
is similar to the one found in tellurite glasses There 
are Cd and Te atoms each of which occupy two non- 
equivalent sites, and O atoms in six non-equivalent sites. 
Both Cd atoms are 6-fold coordinated with Cd-0 bond 
lengths in the 2.197-2.476 A range. The bond angles de- 
viate as much as 31 degrees from the ideal octahedral 
bond angles of 90 and 180 degrees. O atoms always link 
one Te atom with one, two, or three Cd atoms, and the 
average coordination of O atoms is three. 

CdaTeOg has been recently investigated due to its 
semiconducting properties and has been proposed as a 
n-type thermoelectric material^ Its atomic structure^ 
is described as a deformed perovskite-type with pseudo- 
orthorhombic monoclinic space symmetry Y2i/n^ in 
which the B-sites are occupied by Cd and Te, and the 
A-sites by Cd. Adding the oxygen atoms, a total of 6 



non-equivalent positions are occupied. The local environ- 
ments of the non-equivalent atoms are shown in Fig. 
The first coordination shell of Te atoms is composed of 
three pairs of O atoms, each pair at a bond length dis- 
tance of 1.904, 1.924, and 1.948 A. These lengths lie in 
the typical range of paratellurite a-Te02 and tellurite 
glasses i^fi^ Cd atoms occupy two non-equivalent sites. 
One of them is surrounded by three pairs of O atoms at 
bond lengths of 2.211, 2.311, and 2.350 A and bond an- 
gles near 90°. This arrangement is the same as in CdO 
with rocksalt structure, where the Cd-0 bond length at 
room temperature is 2.344 Ai£i The other Cd atom is 
surrounded by four oxygen atoms at 2.237, 2.251, 2.278, 
and 2.297 A. Three of these oxygen atoms lie in nearly 
the same plane as the Cd atom, and the other one lies 
nearly perpendicular. This complex configuration can be 
understood as originated by the electrostatic repulsion of 
the other four O atoms located at 2.596, 2.736, 2.761, and 
3.010 A, which complete a distorted octahedron, but are 
too far apart to be considered as bound. The O atoms oc- 
cupy three non-equivalent crystallographic positions and 
are linked to one Te atom and two or three Cd atoms. 
The two arrangements of one O atom with two Cd atoms 
and one Te atom are nearly planar. 

The greatest difference between CdTeOs and CdaTeOe 
is that tellurium atom is 3-fold oxygen coordinated in the 
first and 6-fold coordinated in the latter. This is due to 
the fact that the content of tellurium in CdTeOs is twice 
that of Cd 3 Te0 6 , if CdTe0 3 is regarded as Cd 2 Te 2 6 . 
Moreover, the Te octahedral coordination is typical in 
compounds where Te has degree of oxidation VI. 

In summary, in both CdTeOs and CdaTeOe, there 
are only Cd-0 and Te-0 bonds. Oxygen atoms always 
bridge one tellurium atom with one, two or three cad- 
mium atoms. In the most typical cases an oxygen atom 
bridges one Te atom with two Cd atoms, with local quasi- 
planar arrangement. 



III. ELECTRONIC STRUCTURE 

We have performed calculations of the electronic struc- 
ture in the framework of the Density Functional The- 
ory (DFT) 52 with the Perdew-ZungeoSi parametriza- 
tion of the Local Density Approximation (LDA). The 
calculations have been performed using the ab-initio 
total-energy program VASP (Vienna ab-initio simulation 
program) i^SifiSiSl The valence electronic eigenfunctions 
are obtained from an expansion in plane waves, while the 
core electrons are replaced by the ultrasoft pseudopoten- 
tials supplied in the VASP package. The cutoff of the 
plane wave expansion depends upon the material com- 
position. We have used the maximal values of the rec- 
ommended cutoffs of the components: 168, 115, and 396 
eV for Cd, Te, and O, respectively. 

To obtain a smooth density of states (DOS), the first 
Brillouin zone has been sampled with k-points gener- 
ated according to the Monkhorst-Pack scheme. We used 
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tion (GGA), both being within 2 % of the experimental 
value 6.481 A£i 

As a complement to the band diagram and the par- 
tial DOS analysis, we use the electron localization func- 
tion (ELF) 58 i 59 i 60 i 61 to obtain information on interatomic 
bonds. The ELF is a measure of the electron localiza- 
tion in space, and helps to visualize the shell structure 
of atoms and solids. The ELF is defined in terms of 
the local electronic density p(r), the excess of kinetic en- 
ergy due to the Pauli exclusion principle T(r), and the 
Thomas- Fermi kinetic energy density T^(r). 



ELF 



T(r) 
T h (r) 



(1) 



In a DFT calculation, T(r) and T/,(r) are evaluated from 
the density and the Kohn-Sham orbitals <fi(r), as 

nr,-i5>««i--§^ o) 

T h (r) = 0.3(37r 2 ) 2 / 3 p(r) 5/3 (3) 

By definition, ELF values range between and 1. An 
ELF value close to 1 at a given point indicates a high de- 
gree of localization of the electrons at that point. For the 
homogeneous electron gas, which is used as a reference 
system in its definition, the ELF is equal to 0.5. 



FIG. 4: Total and partial density of states of CdTe. 



CdTe 



N x N x N grids, with N = 29 for CdTe, N = 33 for 
CdO, N = 10 for Te0 2 , N = 4 for CdTe0 3 , and TV = 6 
for CdsTeOg. To eliminate unphysical oscillations in the 
density of states (DOS) due to the k-space grid, the DOS 
has been broadened by gaussian functions with 0.2 eV 
variance. However, for CdTe a broadening of 0.1 eV was 
used to display the structure of the Cd 4d band. An- 
alyzing the partial DOS we can identify the nature of 
the peaks observed in the DOS, and correlate the fea- 
tures of the total DOS with the crystal structure. The 
partial DOS are obtained projecting the electron wave 
functions onto atomic-like orbitals, and integrating them 
in atom-centered spheres. The radii of these spheres are 
the atomic radii times a factor that equals the sum of the 
volumes of the atomic spheres to the cell volume. Our 
selection of the atomic radii is discussed in the Appendix. 

As a check for the calculations we have performed op- 
timizations of the unit cell volume and the ionic positions 
of CdTe0 3 and CdTe. For CdTe0 3 , the obtained lattice 
parameters were 4 % smaller in the LDA than the exper- 
imental ones, while the ionic displacements are smaller 
than 1 % of the lattice parameters. For CdTe, the theo- 
retical lattice constant is 6.43 A in the LDA, and 6.62 A in 
the framework of the Generalized Gradient Approxima- 



Cadmium telluride is a prototype of II- VI semiconduc- 
tors. It has the zinc-blende structure, where the anions 
arrange in a cubic compact structure, and the cations oc- 
cupy one half of the tetrahedral sites between the anions. 
The band structure of CdTe was determined 30 years 
agoggjgj anc j nag been obtained by ab initio calculation 
many times (see e.g. Refs. 20 22.23 24.64). In Fig. 01 we 
show the DOS. The energies are referred to the valence 
band maximum. The top valence band (VB) and the bot- 
tom conduction band (CB) are mainly composed of Te 5p 
bonding levels and Cd 5p antibonding levels, respectively. 
Below the top VB other VBs are placed, originated from 
the levels Cd 5s, Cd 4d, and Te 5s. It is well documented 
that LDA predictions of the optical bandgap in II- VI 
are significantly reduced compared to the experimental 
values (see Ref . I22I and references therein). We have ob- 
tained a direct bandgap E g = 0.573 eV at the point T, 
which is 1.3 eV smaller than the experimental valued 
Note that the scale of Fig. 0] does not allow one to see 
the CB edge in the total DOS. LDA also predicts the Cd 
4d electron bands to occur about 3 eV higher than the 
experimental values, but still isolated between Te 5s and 
Cd 5s bands. This underestimation of the Cd 4d binding 
energy is due to the unphysical electron self-interaction in 
Cd 4d orbitals, which is not correctly compensated by the 
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FIG. 5: Total and partial density of states of CdO. The zero 
of the energy scale is set at the valence band maximum at 
point L. 



LDA exchange-correlation functionali^Si 2 ^ GGA calcula- 
tions suffer the same shortcoming. However, the topol- 
ogy of the band structure is correctly provided by the 
LDA calculations. The bandgap underestimation has two 
sources. One is the error inherent to the LDA exchange- 
correlation potential, and the other is the discontinuity 
of the exact DFT exchange correlation potential upon 
addition of one electron in insulating solids&S& Calcu- 
lations of the exact DFT exchange-correlation potential 
have shown that nearly 80% of the gap error in Si, GaAs, 
and AlAs is inherent to exact DFT £L However, the calcu- 
lation of GaAs in Ref. |63, which is the most similar case 
to CdTe, placed the Ga 3d electrons in the Ga frozen core. 
When the cation d orbitals are treated self-consistently, 
the unphysical proximity of cation d and anion p bands 
caused by the self-interaction error repels upwards the 
higher valence bands, producing a further decrease of the 
bandgap. 
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FIG. 6: Bandstructure of a-Te02- The fractional coordinates 
of the high symmetry points of the first Brillouin zone are: 
r(0,0,0), M(0.5,0.5,0), Z(0,0,0.5), R(0,0.5,0.5) and X(0.5,0,0). 



B. CdO 

CdO is also a II- VI semiconductor. It has rocksalt 
structure, where the cations are placed in the octahe- 
dral interstitials of the close packed anions. As oxygen 
is more electronegative than tellurium, the Cd-0 bonds 
are more ionic than the Cd-Te ones. The band structure 
that we have obtained is similar to others obtained within 
standard LDA 22 and GGA. 29 Hence, we do not repeat it 
here. It shows an indirect bandgap, with the valence 
band maximum (VBM) at point L, and the conduction 
band minimum (CBM) at point T. LDA and GGA er- 
roneously predict a negative bandgap because it improp- 
erly accounts for the electron self-interaction in Cd 4d 
orbitalsi^ 2 * 2 ^ The unphysical self-interaction in the Cd 4d 
orbitals raises these levels with respect to the experimen- 
tal values. The Cd 4d levels repel the upper O 2p levels, 20 
thus decreasing the bandgap. This is the same effect that 
reduces the bandgap in CdTe and other IIB-VIA com- 
pounds. However, the structure of the VB here obtained 
is similar to that calculated using self-interaction- and 
relaxation-corrected pseudopotentialSjS 2 - where the gap 
problem is corrected and the semiconducting character 
of CdO is restored. 

As can be seen from Fig. the Cd 4d bands are iso- 
lated from the top VBs, which in turn are originated 
from the O 2p levels. Nevertheless, as shown by the par- 
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FIG. 7: Total and partial density of states in a-Te02. 

tial DOS, there is a certain mixing between O 2p and Cd 
4d levels. The mixing of Cd 4d and O 2p levels is the 
origin of the splitting of the Cd 4d band and the peak at 
-3.6 eVi£& At some points in the Brillouin zone, certain 
Cd 4d levels have the same symmetry as the O 2p levels 
and form bonding and antibonding combinations, rais- 
ing the O 2p levels and lowering the Cd 4d ones^ Other 
Cd 4d levels, at the same or at a different k-points, have 
different symmetry than O 2p levels, and originate the 
Cd 4d band at higher energy. The DOS shows this effect 
integrated over the Brillouin zone. 

C. Te0 2 

The band structure and the partial DOS of a-Te02 
are shown in Figs. |5] and [7] We have obtained an in- 
direct bandgap of 2.768 eV, with the VBM located at 
A points k = (±0.36, 0, 0), (0, ±0.36, 0) (in fractional 
units) and the CBM at £ points k = (±0.25, ±0.25, 0), 
(±0.25, ±0.25,0). However, the top VB and the bottom 
CB are very flat and the direct bandgap i?£-£=2.851 eV, 
at k = (±0.2375, ±0.2375, 0) is very close to the indirect 
bandgap. For practical purposes a-Te02 can be consid- 




FIG. 8: Top: Iso-surface of the electron density at a value of 
p = 0.51 A" 3 around Te in a-Te02- Middle: Iso-surface of 
crystal-minus-atomic difference of electron density at a value 
Ap = 0.16 A -3 . Bottom: Iso-surface of the electron localiza- 
tion function at value ELF=0.83. 



ered a direct bandgap material. The flatness of the band 
edges also indicates a high degree of electron localiza- 
tion. The partial DOS shows that contributions from O 
2s and 2p, and Te 5s and 5p orbitals are spread along 
a wide range of energies, thus indicating that partially 
covalent bonds are formed. Also, the Te atom contri- 
butions to the DOS are one order of magnitude smaller 
than O atom contributions, indicating a charge transfer 
from Te atom to O atom. 

Figure [S] shows an iso-surface of the electronic den- 
sity and the ELF around a TeOe cluster inside the crys- 
tal structure of a-Te02- The central Te and the four O 
atoms below it form the above mentioned trigonal bipyra- 
mid. Several features should be noted in this figure: (1) 
directionality of the Te-0 bonds in the tbp; (2) a lone 
pair over the Te atom, opposite to the tbp Te-0 bonds; 
(3) an increment of the charge density around O atoms 
in the plane perpendicular to the O-Te bonds (this fea- 
ture confirms the results of cluster calculations^^); (4) 
there are two regions of high electron localization: one 
over the Te atom, which corresponds to the electrons lone 
pair, and the other around O atoms in the plane perpen- 
dicular to the O-Te bonds. Iso-surfaces for higher ELF 
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(not shown in the figure) concentrate at the Te lone pair. 
A smaller increment of the charge density (Ap ~ 0.05 
A- 3 ) was also found at the Te lone pair region. 



D. CdTe0 3 




D 



FIG. 9: Band structure of CdTeC>3. The fractional co- 
ordinates of the high symmetry points of the first Bril- 
louin zone are: r(0,0,0), A(0.5,0,0), B(0,0.5,0),C(0,0,0.5), and 
D(0.5,0.5,0). 

Figure EO shows the band diagram of CdTcOs. The 
VBM is located at the T point, and the CBM is found at 
k = (-3/8, 1/8, ±3/8) and (3/8, -1/8, ±3/8) (in frac- 
tional units). Therefore, the bandgap is indirect, with 
a value £^=3.028 eV. However, the direct bandgaps 



E g (T) 



3.142 eV and E g {k ) 



3.163 eV are very close 



to the indirect gap. These values are slightly smaller 
than the experimental value 3.9 eV reported for amor- 
phous nearly stoichiometric CdTeC>3 thin films4^ 

The total and partial DOS are shown in Fig. ^| The 
partial DOS reveals that O 2s orbitals combine with Te 5s 
and 5p, while O 2p combine with Te 5s and 5p, and Cd 
4d orbitals. This indicates a complex electronic struc- 
ture, the main features of which are: (1) Cd 4d levels 
concentrate in a single peak in the DOS, with a tail in 
the region of O 2p levels, and no gap in between; (2) O 
2p levels dominate the top VB and also contribute to the 
lower VBs, being determinant in the O-Te bonds ; (3) 
Te 5s and Te 5p contributions are distributed through- 
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FIG. 10: Total and partial density of states in CdTeOa 
(elec/f.u.=electrons per formula unit). The orbital partial 
DOS are plotted for each non-equivalent atom. 



out the VB, and are very depopulated in benefit of the 
O 2p levels; (4) the CB edge is dominated by O 2p lev- 
els, with contributions from Cd 5s and Te 5p levels. We 
conclude that the Cd-0 bonds have a certain covalent 
character, expressed in the mixing of the 4d levels of Cd 
with the 2p levels of O. The ionic character is suggested 
by the octahedral coordination of Cd, and by the spheri- 
cal symmetry of the charge distribution around Cd atoms 
(Fig. Illfl . However, in contrast to CdO, the partial DOS 
of CdTeOa shows that the Cd 4d levels concentrate in 
one main peak, which is slightly broader than in CdO, 
and has a background to the higher valence bands with- 
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FIG. 11: Top: Iso-surface of the electron density at a value 
of p = 0.51 A~ 3 around a OsCdOTeO^ cluster inside the 
CdTeOs crystal. Middle: Iso-surface of crystal-minus-atomic 
difference of electron density at a value Ap = 0.16 A -3 . Bot- 
tom: Iso-surface of the electron localization function at value 
ELF=0.83. 
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FIG. 12: Band structure of Cd 3 Te0 6 . The high symme- 
try points are denoted as in the the orthorhombic system: 
X(0.5,0,0), Y(0,0.5,0), Z(0,0,0.5). 



charge density around O atoms in the plane perpendicu- 
lar to the O-Te bonds (this feature is similar to the case of 
a-Te02) (5) there are two regions of high electron local- 
ization: one over the Te atom, which corresponds to the 
electron lone pair, and another one around O atoms in 
the plane perpendicular to the O-Te bonds. Iso-surfaces 
for higher ELF (not shown in the figure) concentrate at 
the Te lone pair, as in paratellurite. A smaller charge in- 
crement was also found at the Te lone pair region. Note 
that p = 0.51 A -3 is a rather low density (compare with 
Fig. I15fl . For higher values of p, the iso-surfaces are 
nearly atom-centered spheres. 



E. Cd 3 Te0 6 



out gap. This can be understood as an effect of the more 
complex atomic environment in this phase with respect 
to rocksalt-structure CdO. However, the absence of gap 
between the Cd band and the O 2p band may be caused 
by the overall tendency of LDA and GGA to underesti- 
mate the binding energies and the bandwidths. The Cd 
5s levels contribution is present in the conduction band. 

Figure lTTI shows an iso-surface of the electronic density, 
the charge density difference (crystalline minus atomic), 
and the electron localization function (ELF) around 
a cluster OsCd-O-TeC^ inside the crystal structure of 
CdTe03. Several features should be noted in this image: 
(1) spherical symmetry of the electronic cloud around the 
Cd atom; (2) directionality of the Te-0 bonds; (3) a lone 
pair near the Te atom opposite to the Te-0 bonds, i.e., 
over the apex of the pyramid; (4) an increment of the 



FigureH"2lshows the band diagram of CdsTeOe- As can 
be seen, it is an insulator material with a direct bandgap 
Eg = 0.714 eV at the T point. The valence bands are 
very flat, but the lowest CB shows dispersion and is also 
isolated from the upper CBs. From the curvature of the 
conduction band we have estimated an effective mass of 
0.2 free electron masses, which is comparable to that of 
II- VI semiconductors. 

The total and partial DOS of CdsTeOg are shown in 
Fig. EH As in CdTe0 3 , the top VB is built from Cd 4d 
and O 2p atomic orbitals, and the lowest VBs are mainly 
composed of O 2s orbitals with small contributions of Te 
5s and 5p orbitals. The bottom CB has contributions 
from all atoms, the largest of which comes from O 2s and 
Cd 5s. The most significant difference with CdTeOs are 
the triple splitting of the peak associated to O 2s levels 



9 




— i r 

Cd(s) — 
Cd(p) 
Cd(d) 




FIG. 13: Total and partial density of states in CdsTeOe 
(elec/f.u.=electrons per formula unit). The partial DOS are 
plotted for each non-equivalent atom. 
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FIG. 14: Comparative plot of the density of states of CdTe, 
CdO, Te0 2 , CdTeOs and Cd 3 Te0 6 . 



pound studied here. 



between —20 and —15 eV, and the soft CB edge, which 
is associated to the small effective mass. 

We find that the O-Cd bonds are similar to those in 
the other compounds examined. Te atoms are octahe- 
drally coordinated and transfer about 4 electron charges 
to neighbor O atoms. The charge density along the line 
of Te-0 bonds is rather similar to the one in CdTeOs, 
and there is no electron lone pair bond of Tc. The elec- 
tron localization function (not shown here) only points to 
charge localization around O atoms, which is consistent 
with the narrow and isolated O 2s bands. 

Figure^^shows comparatively the DOS in all the com- 



IV. DISCUSSION 

In order to compare our results for the electronic struc- 
ture of Cd-Te-0 compounds, in Tabled we show a sum- 
mary of our calculated bandgaps together with previous 
results. 

The electronic structure of CdO is still not well un- 
derstood due to both theoretical limitations, and exper- 
imental facts that obscure its intrinsic properties. On 
one hand, CdO presents a large defect concentration, 
which produces carrier concentrations as large as 10 20 
cm _3 (2SiS a nearly metallic conductivity, a Fermi level 
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inside the conduction band, 68,69 a blocking of the fun- 
damental absorption edge and a bandgap shrinkage. 68 
To further complicate the interpretation of the optical 
data, two indirect bandgaps have been observed below 
the direct bandgap, which are also highly sensitive to the 
temperature. At the temperature of 100 K the accepted 
values 70 for the two indirect bandgaps are 0.84 and 1.09 
eV, and for the direct bandgap is 2.28 eV. 

Previous band diagram 

calculationa 22 i 25 i 26 i 27 i 28 i 29 i 69 i 71 i 72 and the present 
one provide different interpretations for ordering 
of the indirect bandgaps. Early semi-empirical 
calculations^ 9 *^ and correlated Hartree-Fock self- 
consistent calculations 2 ^* 2 ^ located the VBM at a E 
point, while other semi-empiricalpi LDA) 2 ^ 7 . periodic 
Hartree-FockjS 7 . and screened-exchange LDA 28 calcula- 
tions, as well as this work, locate the VBM at point L. 
In all cases the energy differences between the valence 
band maxima at L and at E are small. Hartree-Fock 
calculations^SSiS 7 . overestimate the bandgaps and 
widths, while LDA 22 i 27 i 28 and GGA 29 underestimate 
them. Moreover, the smallest indirec t ban dgap is 
negative in some DFT calculations (Refs. I22l29l and the 
present one). Although the underestimation of DFT 
bandgaps is known to be caused by the self-interaction 
of the Cd 4d levels and the repulsion from the higher 
levels of the T point, it is not easy to assess the 
reason for the differences between the standard-LDA 
calculations. A common feature of the negative gap 
calculations is the use of pseudopotentials (Refs. |22 
and this work) or supposedly the use of a frozen-core 
electron configuration— The positions of Cd 4d bands 
have also been a point of interest. While the early XPS 
measurement assigned to this band a mean binding en- 
ergy of 12.8 eV— with respect to the (undefined) Fermi 
level, recent measurements have corrected this value to 
10.7 eV, or 9.4 eV with respect to the VBM. 2 ! The best 
approximation to this experimental energy up to date is 
provided by the periodic HF calculation of Ref. |2?]- LDA 
calculations underestimate these binding energies, in 
this study giving 6.6 eV. Our LDA calculation provides 
a Cd 4d band splitting equal to 0.9 eV, which agrees 
with the experimental values of 0.87— and 0.8 eV— 

The electronic structure of a-Te02 provided by our 
calculation is compared in the following with previous 
findings. The band diagram of idealized T3O2 with ru- 
tile structure has been calculated using the empirical 
tight-binding method^ and the scalar-relativistic linear 
muffin-tin LDA method— Both calculations show im- 
portant differences in the DOS, particularly the self- 
consistent one^ 6, where a metallic band structure was ob- 
tained. We performed a test calculation with the struc- 
ture of rutile and we also found a metallic DOS. The 
failure of the rutile model of Te02 is due to the asym- 
metry of the Te-0 bonds in the first coordination shell 
of Te in the paratellurite structure, as has been demon- 
strated in clusters calculations^Sii^*^ and in our results. 
An XPS spectra was reported and interpreted with the 



help of a self-consistent-charge discrete variational Xa 
calculation of TeOg clusters— Our findings of the crys- 
tal electronic structure calculation support this interpre- 
tation and our DOS corresponds fairly well with the dis- 
crete levels of the clusters. The most significant differ- 
ence is the existence of a double peak centered at —11 
eV and —9 eV in our DOS, which is composed of mixed 
Te 5s and O 2p orbitals. This feature corresponds to 
the band C of Ref. H^, which is not split in the calcula- 
tions nor in the observed XPS spectrum. Bcrthcreau et 
alll have reported Hartree-Fock calculations of charged 
clusters TeO^ - , which predict the higher occupied molec- 
ular orbital (HOMO) to be antibonding and composed 
of mixed Te 5s-0 2p orbitals and have suggested that 
this molecular orbital corresponds to the Te lone pair. 
This level can be correlated to the conduction band edge 
predicted by our crystal calculation. The crystal calcula- 
tion is more consistent with less-charged clusters, where 
the HOMO of Berthereaud et al becomes the lowest 
unoccupied molecular orbital, and the HOMO is com- 
posed of non-bonding O 2p levels^ in correspondence 
with our valence band edge. Based on their study of 
charged TeOg~ clusters, Suehara et a& estimated the 
net charge n of the TeOg structural unit to be 6—, and the 
net charge of Te atoms to be 3+, This coincides closely 
with our estimation based on ab-initio atomic radii (see 
Appendix). The optical bandgap provided by our LDA 
calculation near 2.8 eV can be compared with available 
data. Uchida— measured the absorption in the pre-edge 
region, in which an Urbach tail behavior can be observed. 
Although no attempt was made to estimate a direct or 
indirect gap, it is clear from their spectra that the funda- 
mental edge of paratellurite lies in the region 3.6-3.9 eV. 
Mansingh and Kumar— measured the optical properties 
of Te02 thin films prepared by vacuum evaporation and 
found an indirect gap at 3.75 eV and a direct gap at 4.2 
eV related to a strong absorption. Although it is not 
clear if their as-grown films were amorphous or polycrys- 
talline, these values provide an indication of the bandgap 
for ideal crystals. Hence, the LDA calculation underesti- 
mates the paratellurite bandgap by approximately 1 eV. 

The theoretical electronic properties of CdTeOa and 
CdaTeOg can be discussed only with reference to amor- 
phous thin films or oxidized CdTe crystals, as there is no 
data for crystalline CdTeOa and CdsTeOg. The bandgap 
of nearly stoichiometric amorphous CdTeOa is 3.9 eV— 
Hence, LDA underestimates the bandgap in 0.9 eV, 
which is the same trend as in Te02, CdO, and CdTc. 
Based on this trend we predict a bandgap around 1.7 eV 
for CdsTeOg- This material has recently been found to be 
semiconducting and susceptible to n-doping— These ex- 
perimental findings are in agreement with our predicted 
gap, conduction band and effective mass. 

Our final comment is devoted to the spin-orbit cou- 
pling, which is known to have a significant effect on the 
bandgaps of M-Te compounds (M=Zn,Cd,Hg) (see e.g. 
Ref. |2CI) . The spin-orbit splitting of the fundamental 
bandgap of these compounds arises mostly from the Te 
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TABLE I: Calculated and experimental values of the fundamental bandgap (in eV) in all the compounds, and Cd 4d energy in 
CdO relative to the valence band maximum at point L. 



CdTe 












E g ,r 


Method 


E 3 ,r 


Method 


E g ,r 


Method 



1.60 (1.92) 

1.59 

1.44 

0.47 

0.18 



TPA— a 

NLPRH 
LAPW2S 
LAPW-SR2S 
LAPW-R 20 



0.3 
0.8 

0.47 (0.79)" 

0.51 

0.29 



PR££ 

SIRC-PPS 
LMTO-R2 
LMTO-SR2 5 . 
LMTO-RS 



0.80 
1.76 



(0.3) 0.573 



PP,GW2& 



ppo 



CdO 

tp i n d 
Eg,l 


tp i n d 
Eg,2 


jpdir 
Eg 


E{CA 4d) 


Method 


0.84 


1.09 


2.28 


-9.4 


TRffi 
XPS££ 



Eg,L-T 


Eg,T,-T 


Eg,T-T 


E{Cd 4d) 


Method 


0.8 


1.2 


2.38 


-6.1 


APWii 


1.18 


1.12 


2.36 


-5.7 


LCAOS 


1.11 


0.95 




-4.5 


APW 69 


6.07 




6.56 


-12.1 


C-HF 2 " 


10.9 




11.0 


-13.2 


AE-HB 2 ^ 


0.5 


0.4 


0.8 


-14.4 


C-HF 25 


-0.6 






-6.2 


pp22 


1.7 


1.7 


3.4 


-8.2 


SIRC-PPS 


10.41 


10.44 


11.66 


-10.3 


AE-HF^i 


0.39 


0.47 


1.61 


-7.1 


AE-LDAS 





0.2 


2.36 




sx-LDAS 


-0.5 




0.7 


-6.5 


LAPW-GGA 29 


-0.465 


-0.386 


0.712 


-6.6 


ppc 



Te0 2 

E 



Method 



E a 



Method 



CdTeOa 
E a 



CdsTeOe 
Method E g Method 



3.6-3.9 
4.2 (dir) 
3.75 (ind) 




OPAS 

0pA 79<i 

LMTO-SR 46 



1.7 (A-r) 
1.9 (r - r) 

2.768 (A - E) 
2.846 (S - E) 



LCAC^ 5 . 
LCAOi 5 - 

PP C 
PP C 



3.9 



3.028 



OPA 48d 



ppc 



1.7 guess e 



0.714 PP C 



TPA: Two-photon absorption at T=4.2 K. 

OPA: One photon absorption at room temperature. 

TR: Thermoreflectance at T=100 K. 

PP: Self-consistent LDA with pseudopotentials, which include scalar relativistic effects. 

NLPP: Empirical non-local pseudopotentials. 

LAPW: Self-consistent LDA linearized augmented plane waves. 

LMTO: Self-consistent LDA linearized muffin-tin orbital. 

NR: No relativistic, SR: Scalar-relativistic (without spin-orbit coupling), R: Relativistic. 
AE: All electrons calculation. 

SIRC-PP: Self-consistent LDA with self-interaction and relaxation corrected pseudopotentials. 
APW: Non-self-consistent augmented plane waves. 
LCAO: Empirical linear combination of atomic orbitals. 
C-HF: All electrons correlated Hartree-Fock. 
sx-LDA: Screened-exchange LDA. 



"The value in () is the average of the IV — Tg bands split by the 
spin-orbit coupling. 

6 The value in () is the estimated correction due to the spin-orbit 
coupling. 

'Present calculation. 

d In amorphous thin films. 

'Theoretical value plus 1 eV, following the overall trend. 
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5p orbitals. As our calculation does not include this ef- 
fect, our calculated gaps should be regarded as averaged 
values. However, in the compounds other than CdTe the 
partial DOS of Te 5p levels is largely reduced, and the O 
2s and O 2p orbitals give the greatest contributions to the 
DOS. Hence, we expect that the spin-orbit effects on the 
bandgap should be much smaller than in CdTe. However, 
calculations with the spin-orbit coupling included could 
be needed in case that precise measurements of the band 
properties were produced. 
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APPENDIX A: DETERMINATION OF THE 
ATOMIC RADII 




FIG. 15: Ionic pseudo charges along bonds. 



To obtain the partial DOS it is necessary to define 
atomic spheres for the integration of the local DOS. The 
atomic radii cannot be unambiguously defined, being just 
measures of their spatial extension, and depend on the 



environment. Partial atomic charges can be obtained by 

TABLE II: Net atomic charges in CdTeO compounds, deter- 
mined through a partial DOS analysis. 



Cd Te O 



Compound 


R(A) 


q 


R(k) 


q 


R(A-) 




CdTe (C) 


1.48 


-0.46 


1.35 


+0.46 






(A) 


1.29 


+0.46 


1.52 


-0.46 






CdO (I) 


1.09 


+1.6 






1.28 


-1.6 


(A) 


1.20 


+1.1 






1.14 


-1.1 


Te0 2 (I) 






0.80 


+4.4 


1.24 


-2.2 


(A) 






0.9 


+3.1 


1.10 


-1.55 


CdTe0 3 (I) 


1.09 


1.45 


0.80 


+4.42 


1.24 


-1.96 


(C) 


1.48 


-1.05 


1.35 


-0.65 


0.73 


+0.57 


(A) 


1.23 


+0.7 


0.9 


+3.3 


1.10 


-1.35 


Cd 3 Te0 6 (I) 


1.09 


+1.6 


0.70 


+5.6 


1.24 


-1.7 


(A) 


1.23 


+0.97 


0.90 


+3.95 


1.1 


-1.14 


(A) Ab-initio. (C) Covalcnt. 80 


(I) IonicM 





integration of the partial atomic DOS up to the Fermi 
level. These partial charges have no exact meaning, but 
can be used as a reference to compare the chemical bonds 
in the different compounds. 

Usually, tabulated covalent and ionic radii are em- 
ployed for the partial DOS calculations (see Table ITT)l . 
Additionally, we have obtained ab-initio radii analyzing 
the density of charge along the bonds. Figure IT51 shows 
the electronic pseudocharge density along the bond di- 
rections for all types of bonds present in the Cd-Te-0 
system. It can be noticed the difference between the Cd- 
O, Cd-Te and Te-0 bonds. While the Cd-Te and Cd- 
O bonds show clear interatomic regions with minimal 
charge density, the Te-0 bond shows an intermixing of 
the atomic electron bonds. For Cd-0 and Cd-Te bonds, 
we estimated the atomic radius as the distance between 
the nucleus and the minimum of the valence electron den- 
sity along the bond. In the case of Te-0 bonds, the elec- 
tronic density grows continuously from Te to O. Hence 
we set the Te-0 boundary at the point where the second 
derivatives is null. As several distances exist in materi- 
als other than CdTe and CdO, we averaged the obtained 
atomic radii. Table ITT1 shows the atomic charges assum- 
ing these ab-initio radii, as well as ionic and covalent 
radii. As can be seen, with the covalent radii, the atomic 
charges in CdTe and CdTeOa have wrong sign. This is 
rather unexpected in CdTe, where the bond length is 
well approximated by the sum of the covalent radii of Cd 
atom and Te atom. Using the ionic radii, the charges 
in Te02, CdTeOa, and CdsTeOe are in better agreement 
from values obtained from a chemical point of view, but 
the charge of Te atom is extremely high, and O appears 
over-ionized in Te02 . In this sense, the charges estimated 
with the ab-initio radii seem undoubtedly more reason- 
able. The change observed in the partial DOS calculated 
with different atomic radii is fundamentally a scale fac- 
tor. 
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